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Abstract. Algorithms for determining quality/cost /price tradeoffs in saturated markets are considered. 
A product is modeled by d real-valued qualities whose sum determines the unit cost of producing the 
product. This leads to the following optimization problem: given a set of n customers, each of whom 
has certain minimum quality requirements and a maximum price they are willing to pay, design a new 
product and select a price for that product in order to maximize the resulting profit. 

An 0(n log n) time algorithm is given for the case, d = 1, of linear products, and 0(n(log n) d+1 ) 
time approximation algorithms are given for products with any constant number, d, of qualities. To 
achieve the latter result, an 0(nfc d_1 ) bound on the complexity of an arrangement of homothetic simplices 
in M. d is given, where k is the maximum number of simplices that all contain a single points. 



1 Introduction 



Revealed preference theory [TOj is a method of determining a course of business action through the 
review of historical consumer behaviour. In particular, it is a method of inferring an individual's or a 
group's preferences based on their past choices. The marketing mix [9] of a product consists of the 4 
Ps: Product, price, place, and promotion. In the current paper, we present algorithms for optimizing 
the first two of these by using data about consumer's preferences. That is, we show how, given data on 
consumer preferences, to efficiently choose a product and a price for that product in order to maximize 
profit. 

Refer to Figure [TJ A product P = (p, q\, . . . , q c i) is defined by a real- valued price, p, and a number 
of real- valued orthogonal qualities, qi,...,qd- The market for a product is a collection of customers 
C = {C\, . . . , C n }, where Cj = (p,, i, . . . , qi t d)- A customer will purchase the least expensive product 
that meets all their minimum quality requirements and whose price is below their maximum price. 
That is, the customer Cj will consider the product P — (p, q±, . . . , q&) if p < pi and qj > q^j for all 
j G {1, ...,d}. The customer Cj will purchase the product if has minimum price among all available 
products that Cj considers. 

We consider markets that are saturated. That is, for every customer Cj there is an existing 
product that satisfies Cj's requirements and among all products that satisfy C,'s requirements, Cj will 
choose the least expensive product. From the point of view of a manufacturer introducing one or more 
new products, this means that all customers are Pareto optimal, i.e., there are no two customers Cj and 
Cj such that qij~ > qjj. for all k S {1, . . . , d} and pi < pj. This is because, in a saturated market, Cj and 
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Figure 1: A sample market with d = 1 and n = 3. A customer will consider any product that is in their 
upper left quadrant. 



d will purchase the same product, namely the lowest-priced product that satisfies d's (and therefore 
also Cj's) requirements. When modelling a saturated market, there is no need to explicitly consider 
existing products since these can be encoded into the customers themselves. 

As an example, consider a market for computers in which an example customer Ci may be 
looking for a computer with a minimum of 8 GB of RAM, a CPU benchmark score of at least 3000, 
a GPU benchmark score of at least 2000, and be willing to pay at most $1500. In addition, there is 
already a computer on the market which meets these requirements and retails for $1200. Thus, this 
customer would be described by the vector (1200,8,3000,2000). If a manufacturer introduces a new 
product (1199,8,3500,2000) (a computer with 8 GB of RAM, a CPU benchmark score of 3500 and a 
GPU benchmark score of 2000 retailing for $1199) then this customer would select this new product 
over their current choice. 

By appropriately reparameterizing the axes, we can assume that the cost, cost(P), of manufac- 
turing a product P = (p, <7i, • ■ • , 5d) is equal to the sum of its qualities 

d 

COSt(P) = ^ qj . 

i=l 

The profit per unit sold of P is therefore 

ppu(P) = p — cost(P) . 

In this paper we consider algorithms that a manufacturer can use when a choosing new product to 
introduce into an existing saturated market with the goal being to maximize the profit obtained. More 
precisely, given a Pareto-optimal market of customers M — {Ci, . . . , C n }, each have d qualities, the 
PRODUCTDESiGN(d) problem is to find a product P* G M d+1 such that 

profit (P*) = ppu(P*) x \{i : C t purchases P*}\ 

is maximized. To the best of our knowledge, this is the first time a problem like this has been considered 
from an algorithmic perspective. 
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In the remainder of the paper we give an 0(n log n) time algorithm for ProductDesign(I) 
(Section^, and 0(n(\ogn) d+1 ) time approximation schemes for ProductDesign(c?) (Section [3] and 
Section 2]). Section [5] summarizes our results and concludes with directions for future research. 



2 Linear products 

In this section, we consider the simplest case, when a manufacturer wishes to introduce a new product in 
which the quality of a product has only one dimension. Examples of such markets include, for example, 
suppliers to the construction industry in which items (steel I-beams, finished lumber, logs) must have a 
certain minimum length to be used for a particular application. An overly long piece can be cut down 
to size, but using two short pieces instead of one long piece is not an option. 

Throughout this section, since d = 1, we will use the shorthand P = (p, q) for the product 
being designed and qi for q^i. Thus, we have a set of customers M = {(pi, qi), ■ . ■ , (p n , q n )} and we are 
searching for a point P* = (p*,q*) that maximizes 

profit(p*,g*) = (p* -q*)\{i :p* <Pi and q* > q % }\ . 

Our algorithm is an implementation of the plane-sweep paradigm 3]. The correctness of the 
algorithm relies on two lemmas about the structure of the solution space. The first lemma is quite easy: 

Lemma 1. The value (p*,q*) that maximizes profit (p*, q*) is obtained when p* = pi and q* = qj for 
some i, j £ {1, . . . , n}. 

Proof. First, observe the obvious bounds on p* and q*: 

minjpi : i £ {1, . . . , n}} < p* < maxjpi : £ 6 {1, ... , n}} 

and 

min{(?i : £ £ {1, . . . , n}} < q* < max{gi : i £ {1, . . . , n}} . 

Consider the arrangement of lines obtained by drawing a horizontal and vertical line through each 
customer (pi,qi) for i £ {l,...,n}. Within each cell of this arrangement, the function profit(p, q) is 
a linear function of p and q and it is bounded. Therefore, within a particular cell, the function is 
maximized at a vertex. Since each vertex is the intersection of a horizontal and vertical line through a 
pair of customers, the lemma follows. □ 

The following lemma, illustrated in Figure [21 is a little more subtle and illustrates a manufac- 
turer's preference for lower-quality products: 

Lemma 2. Let q' < q and let p be such that < profit(p, q) < profit(p, q'). Then, for any p' < p, 
profit(p', q) < profit (p' , q 1 ). 

Proof. By definition, profit(p, q) = a(p — q) and profit(p, q') = a'(p — q'), where a and a' are the number 
of customers who would consider (p, q) and (p, q'), respectively. These customers are all taken from the 
set M> = {(p u qi) e M : p t > p). 
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Figure 2: profit(p, g) < profit(p, q') implies that profit (p', q) < profit (p', q') for all p' < p. 

Now, consider the customers in the set M' = {(pj, <7i) £ M : p' < pi < p}. By the assumption 
that customers are Pareto optimal, any customer (p,, qi) in M' has qi < q', so all of these customers will 
consider either (p',q') or (p',q) if either one is offered. Therefore, 



Lemma [2] allows us to apply the plane sweep paradigm with a sweep by decreasing price. It tells 
us that, if a product (p, q') gives better profit than the higher-quality product (p, q) at the current price 
p, then it will always give a better profit for the remainder of the sweep. In particular, there will never 
be a reason to consider a product with quality q for the remainder of the algorithm's execution. 

Let the customers be labelled (pi, gi), ■ • • , {p n , In) m decreasing order of pi, so that pi + \ < pi 
for all i € {1, . . . , n — 1}. At any point in the sweep algorithm, there is a current price p, which starts 
at p — oo and takes on the values p%, . . . ,p n , successively, during the execution of the algorithm. At 
all times, the algorithm maintains a list L of qualities q\ > q% > ■ ■ ■ > q^ such that profit(p, q*) > 
profit(p, g|) > • • • > profit(p, q* n ). The quality q\ is the optimal quality for the current price, p. By the 
time the algorithm terminates, the quality of the globally-optimal solution will have appeared as the 
first element in L. 

To complete the description of the algorithm, all that remains is to show how L is updated 
during the processing of a sweep line event. For this, the algorithm uses an auxiliary structure D 
to efficiently identify items in L that need to updated. Consider a consecutive pair of the elements 
q* and q* +1 in L. When q* and q* +1 became adjacent in L, it was at some price p — pt such that 



profit^, q') = (a' + \M'\)(p' - q') 

= a'(p'-q') + \M'\(p'-q') 

>a'(p'-q') + \M'\(p'-q) since q > q' 

= a'(p-q')+a'(p'-p) + \M'\(p , -q) 

> a'(p - q') + a(p' - p) + \M'\(p' - q) since a>a' and (p' - p) < 

> a(p — q) + a(p' — p) + \M'\(p' — q) by assumption 
= a(p'-q) + \M'\(p'-q) 

= profit (p', q) , 



as required. 



□ 
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profit^, q*) > profit(pt, Let <2j and aj+i be the number of customers who would consider (j>t,q*) 

and (pt, respectively. Then, 

profit (p t ,q*) = (pt - 

and 

profit (p t ,qt+i) = (Pt - g*+i)a l+ i 

Now, looking forward in time to a later step in the execution of the algorithm, when p = p t > , with t' > t, 
we find that 

profit (p t >, q*) = (p t > - gf )(dj +t' - t) 

and 

profit (p t >, q* +1 ) = {p t > - q* +1 )(a l+1 + t' - t) . 

We are interested in identifying when the inequality profit (pt>, q*) > profit (pc, q* + i) changes to profit (j>t> , qf) < 
profit(j>4', ql + i). That is, we need to identify all indices i for which, 

(pv - q*)(di +t'-t)< { Pt , - q* +1 )(a l+1 +t'-t) , (1) 

at which point q* should be removed from L. Observe that the values of <Zj, Oj.fi, q*, <Z* + i, and t are all 
fixed at the time q* and become adjacent in L and the only values that change are those of pt> and 
t'. Thus, IJ) defines a halfplane in the plane with axes p t > and t'. 

The auxiliary data structure D used by the algorithm must therefore be able to store halfplanes 
and handle insertions, deletions, and queries of the form "Given a point (x,y) return all halfplanes that 
contain (x,y)." There are several data structures that solve this problem, but the most suitable for the 
current application is the recent dynamic convex hull data structure of Brodal and Jacob [3]0 Their 
data structure allows for the insertion and deletion of halfplanes in O(logrt) time. Given a query point 
(x,y), the data structure is able to, in O(logn) time, find a single halfplane (if one exists) that contains 
{x,y). 

The data structure D is used as follows. When the sweep line is advanced to a new price p t ' , 
the value q t > is appended to L and the halfplane defined by q t > and its predecessor in L is inserted into 
D. Next, the data structure D is repeatedly queried with the point (pt',f). This returns a halfplane 
h (if any exists) that corresponds to a pair of consecutive elements (q*,q*+i) such that profit (p t >, qf) < 
profit (pf, q* + i). The halfplane h is then deleted from D, qi is deleted from L, and a new halfplane 
corresponding to the (now adjacent) elements qi-i and qt+i is inserted into D. This process is repeated 
until querying D with the value (pt',t') returns no result. 

Note that, after all the processing associated with updating the price pt>, the first element, q\, 
in L is the value that maximizes profit (ft., £>{')• Thus, the algorithm need only keep track, throughout 
its execution, of the highest profit obtained from the first element of L, and output this value at the end 
of its execution. This completes the description of the algorithm. 

Theorem 1. There exists an 0(n\ogn) time algorithm for ProductDesign(I). 

Proof. The correctness of the algorithm described above follows from 2 facts: Lemma Q] ensures that the 
optimal solution is of the form (pi,q*) for some i 6 {1, . . . ,n}, and Lemma [2] ensures that the optimal 
solution appears at some point as the first element of the list L. 

This data structure is infamously complicated. It's use in our application can, however, be replaced with a much 
simpler semi-dynamic data structure [6]. So as not to distract from the problem at hand, we defer the discussion on how 
to do this until Section 151 
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Figure 3: Reducing Element-Uniqueness to ProductDesign(1,1). 



The running time of the algorithm can be bounded as follows: Presorting the customers by 
decreasing order of price can be done in O(nlogn) time using any 0{n log n) time sorting algorithm. 
Each sweep line event involves 1 insertion into D plus some number k of deletions, and insertions, and 
k+ 1 queries. Note that each deletion in D corresponds to a deletion in L, and each element of qi, . . . , q n 
is inserted into L at most once. Therefore, the total number of such deletions during the entire execution 
of the algorithm does not exceed n, and each such insertion/deletion pair takes O(logn) time. Since 
there are n events, we conclude that the total running time of the algorithm is 0(n log n) , as claimed. □ 

The following theorem shows that a running time of Q(n\ogn) is inherent in this problem, even 
when considering approximation algorithms. 

Theorem 2. Let M be an instance of ProductDesign(I) and (p*,q*) be a solution that maximizes 
profit (p* , q*). In the algebraic decision tree model of computation, any algorithm that can find a solution 
(p,q) such that 2 • profit (p, q) > profit (p*, q*) has Q(nlogn) running time in the worst-case. 

Proof. We reduce from the integer Element-Uniqueness problem, which has an il(n log n) lower bound 
in the algebraic decision tree model Given an array A = [xi, . . . , x n ] containing n integers, are all 
the elements of A unique? 

We convert A into an instance of ProductDesign(I) in 0(n) time as follows (refer to Figure [3]). 
For each Xi, i G {1, . . . , n} we introduce a customer (pi, qi) with pi — qi + 1/2 and qi = X{. If there exists a 
value x in A that occurs 2 or more times, then the product (x+1/2, x) gives a value profit (a; + 1/2, x) > 1. 
On the other hand, if there is no such x, then 



1. any product (p, q) with p — q > 1/2 can not be sold to any customers and 

2. any product (p, q) with p — q > can be sold to at most 1 customer. 



Therefore, if all the elements of A are unique, then profit (p*, q*) = 1/2, otherwise profit (p*, q*) > 1. 
The result follows. □ 



G 



3 A near-linear approximation algorithm for bidimensional products 

In this section, we consider algorithms for ProductDesign(2), in which products have 2 qualities. As 
a baseline, we first observe that, if we fix the value of (72, then the optimal solution of the form (p, qi, (72) 
can be found using a single application of the algorithm in TheoremQ] Therefore, by successively solving 
the problem for each (72 € {92.1, ■ ■ ■ , 92, n} and taking the best overall solution we obtain an 0(n 2 logn) 
time algorithm for ProductDesign(2). 

More generally, PRODUCTDESiGN((i) can be solved using 0(n d ~ 1 ) applications of Theorem [T] 
resulting in an 0(n d log n) time algorithm. Unfortunately, these are the best results known for d > 2, 
and, as discussed in Section [5j we suspect that an algorithm with running time o{n d ) will be difficult 
to achieve using existing techniques. Therefore, in this section we focus our efforts on obtaining a 
near-linear approximation algorithm. 

Fix some constant e > 0. Given an instance M of ProductDesign^), a point P S R d+1 is a 
(1 - e)-approximate solution for M if profit(P) > (1 — e) profit(P*) for all P* 6 M d+1 . An algorithm is 
a (high probability) Monte-Carlo (1 — e)- approximation algorithm for PRODUCTDESiGN(d) if, given an 
instance M of size n, the algorithm outputs a (1 — e)-approximate solution for M with probability at 
least 1 — n~ c for some constant c > 0. 

Let r — max{ppu(Ci) : i € {1, . . . , n}} and observe that r is the maximum profit per unit that 
can be achieved in this market. Let E = 1/(1 — e) and let I = [log E n\ and observe that t — 0(e _1 log n)H 
For each i £ {0, 1, 2, ... , £}, define the plane Hi = {(p, q\, 52) : P — qi — <l2 = r (l — e) 1 }- The following 
lemma says that a search for an approximate solution can be restricted to be contained in one of the 
planes Hi. 

Lemma 3. For any product P* = (p*,g|,g|), there exists a product P = (p, (71,(72) such that P G Hi 
for some i £ {0, ■■■,£} and profit(P) > (1 - e) profit(P*). 

Proof. There are two cases to consider. If ppu(P*) < r/n then profit(P*) < r, in which case we set 
P = d where ppu(C 8 ) = r, so that P e H a and profit(P) = r > profit(P*) > (1 - e) profit (P*), as 
required. 

Otherwise, r/n < ppu(P*) < r. In this case, consider the plane Hi where i = \\og E {r/ ppu(P*))] . 
Notice, that for any point P € Hi, ppu(P) > (1— e) ppu(P*). More specifically, the orthogonal projection 
P = (p, (71,(72) of P* onto Hi is a product with p < p*, q\ > q\, and qi > q%. Therefore, any customer 
who would consider P* would also consider P, so profit(P) > (1 — e) profit(P*), as required. □ 

Lemma |3] implies that the problem of finding an approximate solution to ProductDesign(2) 
can be reduced to a sequence of problems on the planes Ho, .... Hg. Refer to Figured! Each customer 
Cj considers all products in a quadrant whose corner is Cj. The intersection of this quadrant with 
Hi is a (possibly empty) equilateral triangle Ajj. The customer Cj will consider a product P in Hi 
if and only P is in Ajj. Thus, the problem of solving ProductDesign(2) restricted to the plane Hi 
is the problem of finding a point contained in the largest number of equilateral triangles from the set 
Ai = : j S {l,...,n}}. 

2 This can be seen by taking the limit lim e _ ( Q+ (e/ log(_E)) using one application of L'HopitaFs Rule. 
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Figure 4: The intersection of Hi with customers' quadrants is a set of homothetic equilateral triangles. 

Note that the elements in Aj are homothets (translations and scalings) of an equilateral triangle, 
so they form a collection of pseudodisks and we wish to find the deepest point in this collection of 
pseudodisks. No algorithm with running time o(n 2 ) is known for solving this problem exactly, but 
Aronov and Har-Peled [T] have recently given a Monte-Carlo (1 — e)-approximation algorithm for this 
problem that runs in time 0(e~ 2 n\ogn). By applying this algorithm to each of Aj for i G {1, . . . ,£}, we 
obtain the following result: 

Theorem 3. For any e > 0, there exists an 0(e _3 n(logn) 2 ) time (high-probability) Monte-Carlo (1 — e)- 
approximation algorithm for PRODUCTDESlGN(2). 



4 A near-linear approximation algorithm for constant d 

In this section we extend the algorithm from the previous section to (approximately) solve PRODUCTDESiGN(d) 
for any constant value of d. The algorithm is more or less unchanged, except that the proof requires 
some new results on the combinatorics of arrangements of homothets. 

As before, let r = max{ppu(Ci) : i G {l,...,n}} and let I — |~log E n~|. For each i G 
{0, 1,2,..., £}, define the hyperplane H l = {(p, q 1} . . . , q d ) : p - Y? l= i <li = r ( 1 _ e T)- Tne following- 
lemma has exactly the same proof as Lemma [3] 

Lemma 4. For any product P* — (p* , q\, . . . , q^), there exists a product P = (p, q±, . . . , qd) such that 
P e Hi for some i G {0, ■■-,£} and profit(P) > (1 - e) profit(P*). 

Again, each customer Cj defines a regular simplex Ajj in Hi such that Cj will consider P G Hi 
if and only if P G A^j. In this way, we obtain a set A^ = {A$ i, . . . , A^n} of homothets of a regular 
simplex in M. d and we require an algorithm to find a ((1 — e)-approximation to) the point that is contained 
in the largest number of these simplices. The machinery of Aronov and Har-Peled [T] can be used to help 
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solve this problem, but not before we prove some preliminary results, the first of which is a combinatorial 
geometry result. 

Lemma 5. Let A be a set of n homothets of a regular simplex in W 1 , for d = 0(1), and such that no 
point in M. d is contained in more than k elements of A. Then, the total complexity of the arrangement, 
A(A), of the simplices in A is 0(nk l ). 

Proof. We first consider the simpler case in which the elements of A are translates (without scaling) of 
a regular simplex. Suppose that the total complexity of A(A) is m. Then, by the pigeonhole principle, 
there is some element T in A that is involved in m/n features of A(A), Note that this implies that T 
intersects all the elements of a set A'CA with |A'| = Vt((m/n) 1 ^ d ^ 1 ^), since otherwise there are not 
enough elements in A' to generate m/n features on the surface of A. 

Observe that, since the elements of A' are all unit size and they all intersect T, that they are 
all contained in a ball of radius O(l) centered at the center of T. Furthermore, since each element of 
A' has volume fi(l) this implies that some point must be contained in Q((m/n) 1 ^ d ~ 1 " > ) elements of A'. 
Thus, we obtain the inequality k > Q((m/n) 1 ^ d ~ 1 ' > ), or, equivalently, m < 0{nk d ~ x ), as required. 

Now, for the case where the elements of A are homothets (translations and scalings) of a regular 
simplex, we proceed as follows. Assume, by way of contradiction, that |^4(A)| > rn for some r to be 
defined later. Label the elements of A Ti, . . . , T n in increasing order of size and consider the smallest 
element Ti such that Ti contributes at least r features to A({Ti, . . . , T n }). Such a T,; is guaranteed to 
exist, since otherwise |^4(A)| <rn. 

Now, Ti intersects all the elements in some set A' C {T i+1 , . . . ,T n } with |A'| = ^(r 1 /^- 1 )). 
Shrink each element T' in A' so as to obtain an element T" such that (a) the size of T" is equal to 
the size of T and (b) T" C T'. Call the resulting set of shrunken elements A". Condition (a) and the 
packing argument above imply that there is a point p G R" 1 that is contained in f^r 1 /^ -1 )) elements of 
A". Condition (b) implies that p is contained in ^(r 1 /^ -1 )) elements of A' and hence also A. Therefore, 
we conclude, as before that r < 0(k d ^ 1 ). Thus, for a sufficiently large constant c, setting r = cfc d_1 
yields a contradiction to the assumption that |A(C)| > rn. We conclude that |A(C)| = 0(nfc d_1 ), as 
required. □ 



Remark. The proof of Lemma [5] makes almost no use of the assumption that the elements of A are 
regular simplices other than using the property that their volume is related to their diameter. Thus, a 
version of this lemma holds for collections of fat objects, a result that will probably be of independent 
interest. 



Remark. Lemma [S] is somewhat surprising, since the union of n homothets of a regular tetrahedron 
in, for example, K 3 can easily have complexity f2(n 2 ). This fact makes it impossible to apply the "usual" 
Clarkson-Shor technique [5] that relates the complexity of the first k levels to that of the boundary of 
the union (the 0- level). 

Lemma 6. Let A be a set of n homothets of the regular simplex in M. d such that no point of M. d 
is contained in more than k simplices of A. Then the arrangement A{A) of A can be computed in 
O {nik 11 - 1 + (log n) d )) time. 
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Proof. Computing the arrangement A(A) can be done in the following way. Sort the elements of A by 
decreasing size and construct A{C) incrementally by inserting the elements one by one. When inserting 
an element T, use a data structure (described below) to retrieve the elements of A that intersect T and 
discard the elements that are smaller than T. The proof of Lemma [S] implies that there will be at most 
0(k) such elements. The intersection of the surfaces of these 0(k) elements with the surface of T has 
size 0(k d ~ 1 ) and can be computed in 0(k d ~ 1 ) time using d+ 1 applications of the standard algorithm 
for computing an arrangement of hyperplanes in R d_1 pj [8]. Thus, ignoring the cost of finding the 
elements that intersect T, the overall running time of the algorithm is 0(nfc d_1 ). 

All that remains is to describe a data structure for retrieving the elements that intersect a given 
simplex T 6 A. In the following we describe a data structure that can be constructed in 0(n(logn) d ) 
time and can answer queries in 0(x + (logn) d ) time, where x is the size of the output. This data 
structure will be constructed once and queried n times. The total size of the outputs over all n queries 
will be the 0(\A(C)\) = Oink^ 1 ). 

Since the elements of A are homothets, each element of A can be described concisely by the 
lexicographically smallest of its coordinates and one real number describing its size. In this way, each 
element of A can be represented as a point in M. d+1 and thus A can be represented as a set A' C M. d+1 
of n points . 

Observe that, given a half space h in M. d , the set of elements of A intersected by h is given by 
the set of points of A' contained in some halfspace h! — f(h) in Furthermore, the transformation 

/ preserves parallel halfspaces. That is, if h\,hi Q K d are halfspaces that have the same inner normal 
then f(hi) and /(/12) also have the same inner normal. 

Now, observe that any simplex T G A can be expressed as the intersection of d closed halfspaces, 
all of whose inner normals are taken from the set of N — {n\, . . . , nd+i] C § d_1 of inner normals of the 
facets of a regular simplex. Let f(rii) G § d denote the inner normal of a halfspace f(h), where h C M. d 
is a halfspace whose inner normal is rij. 

The data structure for storing A' is a d + 1 layer range tree [2] in which the ith layer, for 
i G {1, . . . , d + 1}, stores its points ordered by their projection onto f(ni). In this way, the range tree 
can return the set of all simplices in A that intersect T. The size of this range tree is 0(n(logn) d ) and 
it can answer queries in time 0(x + ?i(log n) d ) where x is the size of the query result. Since each simplex 
in A is passed as a query to this data structure exactly once, the total sizes of outputs over all n queries 
is equal to the number of pairs T\, . . . , T% G A such that T\ PI Ti ^ 0. But the number of such pairs is 
certainly a lower bound on |A(A)| so it must be at most 0[nk d ~ ls ). This completes the proof. □ 

Lemma|6]can be used as a subroutine in the algorithm of Aronov and Har-Peled [U Theorem 3.3] , 
to obtain the following Corollary. 

Corollary 1. Let A be a set of n homothets of a regular simplex in R d such that some point p G K d is 
contained in S elements of A. Then there exists an algorithm whose running time is 0(e~ 2d n(\og + 
n(logn) d ) and that, with high probability, returns a point p 1 G M. d contained in at least (1 — e)S elements 
of A. 

As before, an approximate solution to PRODUCTDESiGN(<i) problem reduces to finding deepest 
point in each of the sets Ai, . . . , Ag where Aj is a set of n d-simplices in Hi. By using the algorithm of 
Corollary [T] to do this we obtain the following result: 
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Theorem 4. For any e > 0, there exists an 0(e ( 2<i+1 )n(logn) d + n(logn) d+1 ) time (high-probability) 
Monte-Carlo (1 — e)- approximation algorithm for PRODUCTDESlGN(d). 



5 Conclusions 



We have given an 0(n log n) time exact algorithm for solving ProductDesign(I) and 0(n(log n) d+1 ) 
time approximation algorithms for solving PRODUCTDESiGN(d). The running time of the exact ProductDesign(I) 
algorithm is optimal and no algorithm that produces a (2 — e)-approximation, for any e > 0, can run in 
o(nlogn) time. 

In developing these algorithms, we gave a proof (the proof of Lemma [SJ that shows that an 
arrangement of n fat convex objects in R d has complexity 0(nk d ^ 1 ) where k is the maximum number 
of objects that contain any given point. We expect that this result, and the algorithm for approximate 
depth that arise from it pQ, will find other applications. 

Our exact algorithm for the case d — 2 uses the recent dynamic convex hull data structure of 
Brodal and Jacob [4], which is quite complicated. We observe that the use of this structure can be 
avoided by using the (much simpler) semi-dynamic data structure of Dobkin and Suri [6] . To use this 
structure, we need to specify, each time a point is inserted, the time at which that point will be deleted. 
In our case, Lemma [5] implies that this deletion time can be computed by a binary search on p t , . . . ,p n 
using U). 

An exact near-linear time algorithm for the case d = 2 seems to be out of reach. It appears as 
if this problem requires (at least) a solution to the problem of finding a point contained in the largest 
number of homothets of an equilateral triangle, a problem for which no subquadratic time algorithm 
is known. Is it possible to prove some kind of a lower bound? The related problem of finding the 
point contained in the largest number of unit disks is 3-Sum hard pQ providing some evidence that this 
problem will be difficult to solve in subquadratic time. 

In this paper we considered the case where the problem is parameterized by the number, d, 
of orthogonal qualities that a product may have. Another case to consider is the case in which a 
manufacturer wishes to introduce some number, k, k > 1, of new products into a market. Is this 
problem NP-hard? Does it have a polynomial time approximation algorithm? 
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